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Abstract 

In the context of the mixing dynamical systems we present a derivation of the Gaussian ensembles distributions 
from mixing quantum systems having a classical analog that is mixing. We find that mixing factorization property 
is satisfied for the mixing quantum systems expressed as a factorization of quantum mean values. For the case of 
the kicked rotator and in its fully chaotic regime, the factorization property links decoherence by dephasing with 
Gaussian ensembles in terms of the weak limit, interpreted as a decohered state. Moreover, a discussion about 
the connection between random matrix theory and quantum chaotic systems, based on some attempts made in 
previous works and from the viewpoint of the mixing quantum systems, is presented. 
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1. Introduction 


Gaussian ensemble theory emerged from the study of complex nuclei and long lived resonance states in 
the 1950s by Wigner Eb, and later by Dyson §]. Wigner’s central idea was that for quantum systems with 
many degrees of freedom like a heavy nucleus, one can assume that the Hamiltonian matrix elements in a 
typical basis can be treated as independent Gaussian random numbers. The main prediction of this approach 
is that the statistical distribution of spacings between adjacent energy levels obeys universal distributions which 
define the Gaussian Orthogonal Ensemble, the Gaussian Unitary Ensemble and the Gaussian Symplectic Ensemble, 
if the Hamiltonian is invariant under an orthogonal, unitary or symplectic transformation, respectively. In 1984 
Bohigas, Gianonni and Schimt [U formulated their celebrated statement (briefly named as BGS conjecture) 
concerning quantum chaotic systems: Spectral of time-reversal invariant systems whose classical analogue are K- 
systems show the same statistical properties as predicted by Gaussian Orthogonal Ensembles. Moreover, Gaussian 
ensembles proved to be powerful tool to study statistical properties in many applications fl S S. 01 • 

Ergodic hierarchy (EH) classifies the chaos of dynamical systems according to the decay of correlations 
between subsets of th^hase space for large times. X-systems correspond to the Kolmogorov level of the EH. 
Related to this, in a quantum extension of the EH was proposed, called the quantum ergodic hierarchy 

g EH) . which expresses the decay of correlations between states and observables in the asymptotic limit. In 
Ugl the chaotic behaviors of the Gasati-Prosen model [flit] and the kicked rotator 1 13, IlO were interpreted 
in terms of the quantum mixing level. 

Using the idea of ranking chaos looking at the decay of correlations as in [@,0, we perform two previous 
steps to study the Gaussian ensembles from the quantum mixing level. First, we deduce the mixing factorization 
property which expresses the classical mean value of a product of observables as a product of mean values. 
Second, we obtain the quantum analogue of this property in the classical limit and apply it to deduce the 
Gaussian ensembles. In this way, the contribution of the present paper is to show that Gaussian ensembles are a 
natural consequence of quantum mixing correlations in the classical limit. 
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2. Gaussian Ensembles 


Gaussian ensembles describe how the Hamiltonian matrix elements are distributed in a chaotic quantum 
system when the details of interactions can be neglected. The surprising prediction capability of the GE lies in 
the simplicity of the assumptions. If we have a quantum system having a JV x JV-dimensional Hamiltonian, in 
addition to normalization, the two conditions for the probability density function of the 

Hamiltonian matrix elements which define the Gaussian ensembles are (see, for instance, [Il2l. pp. 73, 74] 
and El p. 62]) 


= P(Hii)P(Hi2)• • • 


( 1 ) 


and 


P(H' ,H]2, .. .,H') = P(Hn,Hi2, • • 


( 2 ) 


where the transformed Hamiltonian H' is obtained from the original one H by an orthogonal, unitary or symplec- 
tic transformation according to corresponding Gaussian ensemble. Eq. © simply represents the invariance of 
the density probability P(Hii,Hi 2 , • • under an orthogonal, unitary or simplectic transformation. Eq. ([T]) 

expresses that in the fully chaotic regime of a classically chaotic quantum system, the details of the interactions 
are not relevant so the Hamiltonian can be replaced by a matrix whose elements are uncorrelated. 


3. Mixing correlations 

In ergodic theory, the decay of correlations is the most important feature for the validity of the statistical 
description because different regions of phase space become statistically independent when they are enough 
separated in time. More precisely, if one has a dynamical system (r,p, S, {TJ) where T is the phase space, 
p : S ^ [0,1] is a normalized measure on S, and {TJtet is ^ semigroup of preserving-measure transformations 
(J is t)q)ically the real numbers), then the EH correlation between two subsets A,B c r separated a time t is 
mathematically expressed as 

C(r,A,B) = p(r,AnB)-M(A)M(B) (3) 

The mixing level of the EH corresponds to the situation when 


lim C(TtA,B) = 0 (4) 


for all A, B c r. Several examples like Sinai billiards. Brownian motion, chaotic maps, belong to the mixing level 
satisfying the eq. Q. The Frobenius-Perron operator P^ associated to the transformation Tj gives the evolution 
of any distribution / (i.e. / : T ^ [0, oo] with 11/11 = 1) by means of 


f(q,p)dqdp 


T-,A 


fiq,p)dqdp 

Ja 


V AC r , V f e J 


(5) 


where (q,p) c T. When P^ has a fixed point /*, i.e. P^/* = /*, the following relevant property of mixing systems 
can be deduced. 


Lemma 3.1. (Factorization property) Let /* be a normalized distribution which is a fixed point of the Frobenius- 
Perron operator P^. If 1^^, 1^^, • • •, 1a„ • T ^ K are the n characteristic functions of n subsets A^,... ,A„ C r then 

' ( C 

/*(q>p)iAi(q,p) • ■ • iA„(q,p)dqdp = 

Jr vJr 

Lemma [3TT] implies that the average of a product weighted by a distribution /*(q,p) (that is a fixed point of PJ 
can be factorized in the corresponding product of the averages weighted by the same /*(q,p). The “factorization 
property" of Eq. ® is essential in order to obtain the Gaussian ensembles, we explore its consequences in the 
context of quantum mixing correlations. 


/*(q,p)l^j(q,p)dqdp 


fAASq,P)dqdp 


( 6 ) 
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4. Quantum mixing correlations 

A quantum counterpart of mixing correlation of Eq. ([4]l was derived in [@] . More precisely, in the quantum 
version of Eq. Q we have a decay correlation between states and observables rather than between subsets of 
phase space given by 


c(p(t),d) = {d>p(,)-(o>p^ 


(7) 


where the role played in ® by the subsets A, B is now played by the state p(t) and the observable O, with p(t) 
being any quantum state p at time t. The state p* is the weak limit of p given by the quantum mixing level of 
the quantum version of the ergodic hierarchy, i.e. the quantum ergodic hierarchy (QEH) 

lim C(p(t),d) = 0 for all observable 6 rg') 


We can see the similarity between the mixing correlation and its quantum version (Eqns. Q and ([8])), that is, 
one can obtain one correlation from the other simply exchanging C[T^A,B) by C(p(t),d) and vice versa. Eq. [8] 
describes the relaxation of any quantum state p with a weak limit p* where the relaxation is understood in the 
sense of the quantum mean values, i.e. the decoherence of observables We show that the weak limit 

p* is the quantum analogue of the distribution/* of Lemma l3.ll This is the content of the following result. 

- -it — 

Lemma 4.1. The state p* is a fixed point of the evolution operator = e » where H is the Hamiltonian of the 
quantum system, i.e. U^p^,UI = p*. 


In order to establish a quantum version of Lemma 13.11 we recall some properties of the Weyl S 5 mibol and the 
Wigner function. If A is an operator then its Weyl S 5 mibol 14^ is a distribution function over phase space defined 

byEli3] 


WA<iq,p) 


J K 



. pA 

e ' » dA 


(9) 


In particular, if lg{q,p) is the characteristic function of a subset B of T we will use the Weyl symbol of fig, with 

= V(q,p)eR^ (10) 

The Wigner function is defined by means of the Weyl symbol as 


WA(.q,P) = ^WAiq,p) 


( 11 ) 


A relevant property of the Wigner function is that it allows to express any quantum mean value as an integral in 
phase space [jlbh . in the form 


(6>p 


dqdp W^(q,p)Wo(q,p) 


( 12 ) 


For the Weyl symbol of a product of operators, it can be shown that the following expansion is fulfilled 


^AsCq-P) = WA(.q,p)WBiq,p) + 0(ft) 


(13) 


An important property that we will use can be deduced by the definition of Weyl s)mibol in the classical limit 
h->0. 

Lemma 4.2. Let Vl^(q,p) be the Weyl symbol of an operator A Then in the classical limit of h 0 the Weyl 
symbol ofWi—t) = tJ^AU^ is W 4 (q(t),p(t)), where (q(t),p(f)) = (T^q, T^p) and is the classical evolution given 
by Hamilton equations. That is, 

^t/;At/,(<3’P) = ^A(q(0,p(t)) V (q,p) G , V t G R (14) 
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For quantum mixing correlations the following property in phase space is a consequence of Lemmas 14.11 and 


Lemma 4.3. The Wigner distribution Wp^(^q,p') is a fixed point of the Frobenius-Perron operator associated with 

the classical evolution given by Hamiltonian equations. 

Now joining the previous Lemmas [3.1114.11 [4.21 and 14.31 we show a quantum analogue of the factorization 
property (i.e. Eq. ®). This is one of the main results of the present contribution. 

Theorem 4.4. (Quantum factorization property) Assume one has a mixing quantum system, i.e., the correlation 
C{p(^t), 6) of p{t) with any observable O vanishes for t —> oo. Then, for a set of observables O^,... ,0^ when ft ^ 0 
one has 

(0i---4)p. = (6i>p.---(4)p. (15) 

Proof. In principle, the Wigner property of Eq. (fT^ applied to the product Oj ■ • • and p* gives us 

r 


{OvOn)p, 


Wp.(q,P)Woi-o„(q,P) dqdp 


Applying several times Eq. (fl^ on Oi ■ • • 0„ we have 

Woi-6„(q,p) = Wo^(.q,p) ■ ■ ■ Wpjq.p) + o(ft) 
From Eqns. (flbl) . (fTTl) and since Wp^{q,p)dqdp = 1 it follows that 


(Oi---On)p. = lVp.(q,P)Woj(q,p)---WoJq,p) dqdp + 0(ft) 
Then, in the classical limit ft ^ 0 we can neglect terms of order 0(ft) so (flSl) becomes 


(Oi---6„)p = 


Wp.(q>P)Woi(q>P)---Wo„(q-P) dqdp when ft^O 


(16) 


(17) 


(18) 


(19) 


We can expand Wq^,..., as linear combinations of characteristic functions. 

That is, Wofq,p) = aylc^.iq,p),... ,Wo„(q>P) = Hi ^hen we have 

(Oi---6„)p. = iv-p.(q,p) ^aijic„(q>p)---^«n(ic„,(q,p) dqdp 

V. ;■ ( J 




Ctnl 


Wp.(q,p)ir, (q>p) • • • ic„,(q>p)‘iqdp 


Now since Wp (q,p) is a fixed point of P^, as shown in Lemma R~T1 then we can apply the factorization property 
of Lemma [3(1] to the integral in the right hand of Eq. (1201) : 


Wp, lc„ (q, p) • • • lc„, (q, p)dqdp = 


Wp,lci^(q,p)dqdp- 


Wp lc„,(q,p)dqdp 


( 20 ) 


This implies that 

Hj «ij ■■■Hi ^nl Jh 2 Wp,(q,p)lcy(q,p) ■ • • lc„,(q,p)dqdp 
= Jh2 WpS^,P)Hj aijlc„(q>P)‘lqdp^p,(q,p)Hi anilc„,(q,P)dqdp 
= Jh 2 Wp,(q,p)Woq(q,p)dqdp ■ ■ ■ Wp^iq,pW 6 S^’P^dqdp = 
(Ol)p.---(On)p. 

which ends the proof. 

Theorem 14.41 expresses the quantum version of mixing correlations in the classical limit ft —» 0. 


( 21 ) 

□ 
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5. Gaussian ensembles by means of mixing quantum systems 


The manifestation of chaotic aspects in quantum systems is possible within characteristic timescales t 


< 


(with T (X in the regular case being a proportional to the phase space dimension, and t oc —logh in the 
chaotic case). In these timescales, the semiclassical and quantum descriptions overlap with the particularity that 
in the logarithmic timescale — log ft the statistical predictions of the Gaussian ensembles are displayed [11 


13 


18 . 1^ . Moreover, within the logarithmic timescale it is expected that the states contain statistical properties as 


the randomness and invariance conditions [Eqs.([T]) and @], and expressed in terms of quantum correlations. 
This motivates the following connection between Gaussian ensembles and mixing quantum systems. 

As we have shown in Sec. 4, the quantum correlations of mixing quantum systems are contained in the 
weak limit p* which is representative of the quantum system in the asymptotic limit t ^ oo. And since the 
logarithmic timescale imposes that t < — log ft then the asymptotic limit can be guaranteed in the classical limit 
for ft vanishingly small. 

The following lemma constitutes a useful tool in order to deduce the Gaussian ensembles within the mixing 
quantum formalism. 


Lemma 5.1. Assume one has a quantum system S subject to a Hamiltonian H with Hamiltonian matrix ele¬ 
ments having a density probability function P{Hii,Hi 2 , ■ ■ ■ Let • j^nnC-^nn) be the 

marginals of P{Hii,Hi 2 , ■ ■ ■ ,Hnn') with respect to the variables Hii,Hi 2 , ■ ■ ■ ,Hnn- Then for each set of values 
■ ■ ■ ,HnnI PiiiHii), PiiiHii), ■ ■ ■,PNNiHNN) and for dHjj, • • •, dH^^^ sufficiently small, there ex¬ 
ist projectors 'ft{Hii,Hi 2 , ■ ■ ■ ,Hnn\ ^ 12 (-^ 11 j^ 12 j ■ • • ■ V ^NNCfflljffl 2 j • • • jffNwX 

and a weak limit pffHii,Hi 2 , ■ ■ ■ ^ach that 


^ i,j = l,...,N (22) 

and 

(fc(Hii,Hi2, . . . ,Hnn'))— Pi^lliHlh ■ ■ ■ (23) 

Moreover, in the classical limit h—>Q, the projector n(Hii,Hi 2 , ■ ■ • ^an be expressed in terms of 

^1i(-^11j-^12j • • • j-^NnX ^12(.^ll>bii2’ ■ ■ ■ ■ ■ ■ > ^Nn(-^11j-^12j ■ • • aS 


^ — ^11^12 ■■■ ^NN- (24) 

Two remarks can be made regarding this lemma. First, the product PiffH^jjdH^j gives the probability that the 
ij-th Hamiltonian matrix element belongs to the interval [Hi^Hij + dH^jj, and a similar statement follows for 
the product P[Hii,Hi 2 , ■ ■ ■ ,H^jv)dHiidHi 2 ■ • • dH^^^ which is the joint probability of the former. Besides, due to 
previous remark and since ft, ftii7ti2, • • •, Unn projectors, then Eqs. (| 22]) and (|23]) express a sort of Born rule 
[f2^ performed by means of weak limit states. Lemma 5.1 allows one to obtain the randomness and invariance 
conditions that define Gaussian ensembles. This is the content of the following theorem. 

Theorem 5.2. (Gaussian ensembles distributions from mixing quantum systems) 

(i) Assuming that S is a mixing quantum system, then in the classical limit ft ^ 0 one obtains the randomness 

CONDITION 


P(Hn,Hi2,... = Pii(Hii)Pi2(Hi2) • • • Pnn^H^n) 

(ii) Let us consider the transformed variables H'-^^,H'-^^, .. • corresponding to Hii,Hi2, ■ ■ • jffNiv through the 
change of variables H' = UHU^ where U stands for the transpose, complex transpose, or dual of U if U 
is orthogonal, unitary, or simplectic, respectively. Let P(H^j,H^ 2 j • • • be the transformed probability 

density function of P { Hii , Hi2 , ■ ■ • .^nn)- One obtains the invariance condition 

P(H'„H' 2 , ... ,H;^) = P(Hn,Hi 2 ,... 
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Proof. (0 If one applies the quantum factorization property (Theorem I4.4D and the Lemma ITT] to the projec¬ 
tors in the classical limit ft ^ 0 one obtains 

(^11^12 ■ ■ ■ ^Nn)p, = (^ll)p,(^12)p, ■ ■ ■ (^Nn)p, (25) 

where for the sake of simplicity we have omitted the explicit dependence on in all the 

expressions. Now using the Eqs. (|22T) . (1231) and ( [24l) one can recast (1251) as 

... ,H^w)dHiidHi2 ■ • ■ = P^,[H^,)dHuPuiH^2)dH,2 ■ ■ ■ PNN(.Hj^f^)dHf^f., 

Then, since dHiidHi 2 ■ ■ ■ dHi.ji.j are arbitrary small then it follows the desired result. 

(ii) By the Lemma (13.ID there exist projectors 7i{Hii,Hi2, ■ ■ ■ ,Hnn)> ■ ■ ■ ’^nn'^ weak limit 

states p*(Hii,Hi 2 ,...,Hjvjv), such that 

(^)p, = PiHiiiHii! ■ ■ ■ >HNN)dHiidHi2 ■ • • dHf.ij^ 

{Ti')p- = P{H[^,H[^, .. .,H'^^)dH[^dH[^ ■ ■ ■ dH'^^ (26) 

where again we have omitted the explicit dependence on Hii,Hi 2 , ■ ■ -^Hnn- Since p' and ft' refer to the 
transformed density probability then it must be satisfied that 

7i' = UnU^ , p'=i/p,U^ (27) 

From Eqs. (l26l) and (|271 ) one obtains 

P(H[^,H[^, .. .,H'^^)dH[^dH[^ ■ ■ ■ dH'^^ = {7f')p- = {U7iU-^)op,D^ = TriUp^UWum = 

= P(Hn,Hi2,... ,H^w)dHndHi2 ■ • ■ dHjvw (28) 

where Tr(...) stands for the trace operation. Since the volume element dHiidHi 2 ■ ■ • dffjvjv is invariant un¬ 
der the transformation, i.e. dH'jdH'^ • ■ • “ dHiidHi 2 ■ ■ ■ then the desired result is obtained 

straightforwardly. 

□ 


6. Physical relevance 

6.1. Kicked rotator 

We illustrate the role played by the Gaussian ensembles in mixing quantum s yste ms with an emblematic 
example of the literature: the kicked rotator 13]. The Hamiltonian is given by fll3. p. 9] 


H = + XcosO 


S 


5(t — uit) 


which describes the free rotation of a pendulum with angular momentum L, periodically kicked by a gravitational 
potential of strength A. The moment of inertia I is normalized to one, and t is the kicking period. We focus on 
the fully chaotic regime, that corresponds to A > 5 G1 pp. 10, 11]. We show that the kicked rotator behaves 
like a quantum mixing system in this regime. Let p be the state of the system at t = 0 

P ='^Pkk\o-k){ak\+'^'^Pki\o-k){ai\ , pfc; = (a^lpla,) Vfc,l (29) 

k ki^l I 

Here {|a^.)} is the Floquet eigenbasis [H^ . with eigenvalues where are the so called Floquetphases. 

For an observable O, after M kicks one has 

(0)p(M.) = Tr(p(MT)(5) = 2p,,0,fc + 22p,,0fc,e-'“(^^-^'^ , O,, = {a,|(5K> Vfc,l (30) 

k ki=l Uk 
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It is shown that for A > 5 the quadratic mean value of the momentum, (L^), exhibits exponential localization 
having a characteristic macroscopic width 1^. Moreover, if M » the phases in the factors in (|30]) 

oscillate rapidly in such a way that only the diagonal terms survive; thus 

{0)p(^M^)^'^PkkOkk forM»Zj (31) 

k 


Note that, if we define the diagonal part of p as p* = {a^\, then we have 

{0)p^ ='^kPkkOkk (32) 

Therefore we deduce that 

Jin^(0>p(MT) = ^ PkkOkk = {0)p, (33) 

Recalling Eqs. (0 and ([8]), it follows that for A > 5 the kicked rotator behaves as a mixing quantum system. 
Even more, p* is the weak limit which is also a mixture of pure Eloquet eigenstates {|a;t)} ^^d then it can be 
interpreted as a decohered state, diagonal in the Eloquet basis, with a decoherence time tp = Nt ~ tIj. 

Summing up, for A > 5 the kicked rotator is a mixing quantum system (with decoherence in the Eloquet basis 
and induced by dephasing) and therefore, the validity of the application of the Gaussian ensembles is justified 
due to Theorem (I5.2D in the classical limit ft ^ 0. 


6.2. Some standard approaches and the mixing quantum formalism 

Here we provide a discussion of the connection between Random matrix theory (RMT) and quantum chaotic 
systems, based on some attempts made in previous works (see, for instance, [1^. |^. 1^ |2^ ) and from the point 
of view of mixing quantum systems. Due to the vast body of work on the subject ( [l25L \2&L 22, 28, ^ ^ 31 


32LI33I1 . among others) and since our starting point are the mixing d 5 mamical systems, we restrict the discussion 


to quantum systems that are chaotic in their classical limit. 

Beyond the success of RMT in the prediction of statistical properties in several phenomena and its consoli¬ 
dation as a specific discipline, there exist questions laying on the foundations of RMT that still remain open or 
partially answered. Below we quote some of these issues and discuss them from the point of view of the mixing 
quantum formalism. 


From the point of view of the quantum mechanics, the redundant information contained in the exponentially 
large number of relevant periodic orbits conceals possible connections between quantum chaotic systems and 
RMT ^fj. 

In the mixing quantum formalism the only relevant information about the system is contained in the 
correlations between the observables and the weak limit, as shown in Eqs. 0 and 0. As we mentioned 
at the beginning of the Section 4, the weak limit process is a t 5 q)e of decoherence of the observables 
[0, 14 , rial . Thus, the redundant information is suppressed by the cancellation of the quantum mixing 


correlations in the asymptotic limit. 


• Although there have been several ways to deduce the BGS conjecture, for instance the non-linear a-model 
1^7 and the semiclassical trace formula ^1^1, all the attempts are based on some kind of semiclassical ap¬ 
proximation 1^7 . 

One of the advantages of the mixing quantum approach is that it allows to deduce the Gaussian ensembles 
distributions as a consequence of the quantum factorization property in terms of operators and states 
(Theorem [50, i.e. in the language of the quantum mechanical operators. However, the classical limit has 
to be considered in order to apply the quantum factorization property. 


• In Gaussian ensembles the behavior is studied along the energy axis rather than the time axis, while the 
thermodynamics systems evolve along the time axis. Furthermore, since there is no way of describing mathe¬ 
matically the transition from one level to the next then there is no analog of the time arrow of thermodynamics 
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My. From Eqs. (|22T) and (|23]) it can be seen that the evolution is involved (in the asymptotic limit) since 
the joint density probability and its marginals 

expressed in terms of traces of projectors in the weak limit, which also allows to give a probabilistic inter¬ 
pretation according to the Born rule. In addition, one can say that in a mixing quantum system the time 
arrow is due to quantum mixing correlations which gives place to an irreversible dynamics expressed by 
the mixture character of the weak limit, as shown in Eqs. (|32]) and (|33]) . 

• Since the trajectories of a dynamical system are system-specific, the role of the ensemble theory of statistical 
mechanics rules out. Instead, ensembles of different Hamiltonians are able to mimic the statistical behavior of 
a dynamical system. 

In the mixing quantum formalism the key point is to consider the statistical description given by the Er- 
godic Hierarchy in terms of correlations between subsets in phase space instead of using trajectories. In 
particular, for the mixing level one has that any two subsets separated enough in time have a null corre¬ 
lation and can be interpreted as statistically independent events. Moreover, this statistical independence 
property can be generalized for subsets in a sequence, the factorization property (Lemma [3jTT) . that char¬ 
acterizes randomness between subsets. In turn, in the classical limit the quantum factorization property 
(Theorem 14.4D allows to express the randomness in terms of the factorization of mean values in the weak 
limit. 

In order to illustrate how the Gaussian ensembles are deduced from the standard approaches in RMT, below 
we provide a schematic picture showing some of them along with the mixing quantum formalism. 


Wigner-Dyson original approach 

universality in local fluctuations of quantum spectra => randomness and invariance Gaussian ensembles 

semiclassical trace formula 

Gutzwiller trace formula => n — point correlation functions in the classical limit => Gaussian ensembles 

mixing quantum formalism 

factorization property quantum factorization property in the classical limit => Gaussian ensembles 


7. Conclusions 

We have proposed a novel way to deduce the Gaussian ensembles within the quantum mixing level of the 
quantum ergodic hierarchy. The relevance of our main contribution. Theorems 14.41 and 15.21 lies in the following 
remarks: 

• In the classical limit the randomness condition of Gaussian ensembles results as a consequence of the 
quantum mixing correlations. 

• The probability density function for the Hamiltonian matrix elements can be computed in terms of the 
mean value of a projector in a weak limit (Eq. (1231) ). In addition, this can be considered as a kind of 
analog of the Born rule. 

• Eor the kicked rotator case we show that mixing quantum formalism links decoherence, in the Floquet 
basis and induced by dephasing, with Gaussian ensembles in terms of the weak limit which also can be 
interpreted as a decohered state. Moreover, starting with a pure state the mixture character of its weak 
limit expresses the irreversible dynamics of the mixing quantum systems, as shown in Eqs. (l32l) and (l33l) . 
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• Going further, from Theorems 14.41 and 15.21 we could rephrase the statement of the Bohigas-Giannoni- 
Schmit conjecture [0] for the family of mixing quantum systems as: Hamiltonian matrix elements of mixing 
quantum systems show, in the classical limit, the same probability density function as predicted by Gaussian 
ensembles. 

Summarizing, we conclude that the “imitation" of statistical properties of quantum systems having a mixing 
(and therefore chaotic) classical limit arises as a consequence of the quantum factorization property within the 
mixing quantum formalism. 
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Appendix A. Proof of Lemma 3.1 

Proof. Let us writeas a linear combination of characteristic functions in the form/*(q,p) = ajlc.(q,p) 
with Q n Cj = 0 if i j and Q = F. Then j^f^,{q,p)dqdp = = 1. For Ai,A 2 Q F, from definition 

dUl we can write 


/^(TfAi n A 2 ) — C{T^Ai,A 2 ) + /4(Ai)p(A2) 
9 


(A.1) 


Let us compute the following expression: 


Hi Hj n (Aj n Cj nA2)) = Xi;«; X;«; [c( 7 ’tCi,Ai n Cj nA2) + /t(Q)p(Cj nAj nA2)] 

= Hi «i Xi; a;C’(rtQ,Ai n Cj nA2) + Xi “iM(Q)Xij “jmCCj nAj nA2) 

= Si “i Sj ajC(rtQ,Ai n Cj nA2) + X; «; Jr 
= Si «i Sj ajC(r^Q,Ai n Cj nA2) + J^X; «;ic/q>p)iAi(q-p)iA2(q.p)'^<j‘ip 

= Si aiXj “jC(rtQ,Ai nCj nA2) +Jj,/*(q,p)l2ij(q,p)lA2(q,P)dq<ip 

Also we have that 


(A.2) 


Si «i Sj ajM( 7 ’tQ nAj n Cj nA2) = Xi «i S; cLji^iJtCi n rtCT.^Aj n Cj nA2) 

= Si “i Sj “jC(Pt(Q n r_tAi), Cj nA2) + Xi “iPCQ n r.tAJXj ajKCj nA2) 

= Si “i Sj ajCCTtCQ n r_tAi), Cj n A 2 ) + Xi «i Jr lQ(q>P)lr_,Ai(q,P)<iqdpXj «j Jr lcj(q.P)lA2(q,P)<iqdp 

= Si aiXjajCCrtCCi nr_tAi),Cj nA2) +J^/*(q,p)lr_,Ai(q.P)dqdp Jr/*(q,P)lA2(q>P)<iq<ip 

= Si «i Sj ajCCTtCQ n r_tAi), Cj nA2) + J^, /*(q,p)dqdp J^/*(q,p)lji^(q,p)dqdp 


Now by the definition of the Frobenius-Perron operator and since /* is a fixed point of P^, we have 


/*(q,p)dqdp : 


f-tAi 


Ptf*iq,P)dqdp ■- 


fM,P)dqdp = 


/*(q,p)iAi(q,p)<iq<ip 


(A.3) 


(A.4) 


Then using ( IA.4D we can recast HA.3D as 

Si «i Sj ajPCT’tQ nAj n Cj nA2) 

= Si «i Sj^jCCTtCQ n r_tAi),Cj nA2) + J^/*(q,p)ljij(q,p)dqdp J^/*(q,p)lji^(q,p)dqdp 


(A.5) 


In the limit t —> 00, the correlations C(rtCi,Ai n Cj nA2) and C{T^{Ci n P.^Aj), Cj nAj) become vanishingly 
small due to Eq. Q since TtCCj n T-tAj) c p^Cj and Cj nA2 are sufficiently separated in time for large times. 
Therefore, from Eqs. (IA.2D and (IA.5D we have 


lim ajpij^Ci nAj n Cj n A 2 ) 


/*(q,p)l, 4 j(q,p)lA^(q,p)dqdp = 


/*(q,p)lAj(q,p)dqdp 


f*iq,P)'i-AXq,P)dqdp 


If we have n characteristic functions 1 ^ 5 ^, 1^^,..., 1 ^ 5 ^ we can apply the last equality n — 1 times so that we prove 
the desired result. □ 


Appendix B. Proof of Lemma 4.1 

Proof. Let s be a real number and let us consider the evolved operator UjOtJ^ for a given operator O. From Eq. 
(H we have 


iim(u;'m),„-(i;; 6 [/,}p,=o (B.i) 

and applying trace properties we can rewrite it as 

Inn - iO)u,p,u] = 0 (B.2) 
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Since 


Im (d}p(,+,) = lim = (6)^^ (B.3) 

then it follows that {6)i; p f/t = (6}p^ for all observable 6, which means that 

U,p,Uj=p, VseR (B.4) 

□ 


Appendix C. Proof of Lemma 4.2 


Proof. From the definition of the Weyl symbol, Eq. one has 


WAiq,p) 


(q + A|A|q — A)e^'^ dA 
Jr 


(C.l) 


Then it follows that 


W^iT.qJ.p) 


r 

{T,q + A\A\T,q - A)e^^~dA 

JR 


Now we make the change of variables A —> A = T_^A, then 

A = TtA and dA=|TjdA 


(C.2) 


(C.3) 


being \T^ \ the Jacobian determinant of restricted to the coordinates q. Moreover, since the Liouville classical 
evolution preserves the volume of phase space we can assume that ITJ = 1. Then, using (IC.3I) and that \T^ \ = 1 
we can recast (IC.2I) as 


WA{T,q,T,p) 


r 

{T^q + T,A\A\T,q - A)e^‘~i~dA 

Jr 


It is clear that 


{T,q + T,A\ = (T,(q + A)| = (q + A\U\j) 
\T,q - T,A) = \Uq - A)) = f/(t)|q - A) 


(C.4) 


(C.5) 


and also 


‘"h = TpTa _ eA — meZ <;=^ p{t)A[t) — pA = mh/2 (C.6) 

Since the quantum phase space is grained due to the Indetermination Principle by cells of volume h/2, then in 
the classical limit for h vanishingly small the condition of (IC.6I) is satisfied. Therefore, replacing (IC.5D and (IC.6D 
in gC.4D we obtain the desired result. □ 


Appendix D. Proof of Lemma 4.3 


Proof. By applying the definition of Frobenius-Perron operator, Eq. ([S]), to the Wigner function VFp/q,p), using 
Lemma RTT] and Eq. (fT^ . we have 


' 

^’tWp,(q>P)dqdp 

Ja 


Wpfq,p)dqdp 

Jr^fA 


Wpfq,p)lT_,Aiq,P)dqdp 

Jr2 


(D.l) 
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Now let ftj^he the operator such that = l/i(q,p). By applying Eq. ( fT^ to it follows that = 

l^(Tfq, Tf-p) = Then using this in Eq. (ID.II) we have 

XiTtWp^(q,p)dqdp = Wp^iq,pWot;^(j^liq,p)dqdp 
= Tr(p*[/J4i/f) = Tr(&tP*LiJ4) = Tr(p*4) = Wp^(.q,PWj^(.q,p)dqdp = W^Jq,p)lA(q,p)dqdp 

= Ia^pM’P^^^I^P ®- 2 ) 

Then, since A c r is arbitrary and given that Pt^p,Cq>p) and Wp^(q,p') are non negative, it follows that 
PtWp^(q,p) = Wp^[q,p) almost everywhere on R^. Nevertheless, since we only use PtWp^[q,p) and W^^(q,p) by 
means of integrals we can consider without loss of generality that PtWp^[q,p) = Wp^(q,p) for all (q,p) e R^. 
This completes the proof. □ 


Appendix E. Proof of Lemma 5.1 

Proof. Since P(Hii,Hi 2 , • •. ,Hjvjv), PuCHu), ^ 12 (^ 12 ), • • •, Pmn(.Hnn) are positive numbers then one can consider 
dHji, dHi 2 , • • •, dHjvjv > 0 sufficiently small such that 

0 < Pi/Hy)dHy <1/2 V i, j = 1,..., JV 


0 < P(Hii,Hi 2, ... ,Hjvjv)dHiidHi2 • ■ • dH^^, < 1/2 
Let a, p be positive numbers such that 

a < min{P(Hj2 ,Hj 2, • • • jTfjviv)?^ii(Tfii)>T’i2(ff22), • • • jPNNi^NN')} 

max {P(Hii, Hi2, PuiHul ^12(^12), ■ • ■, Pnn^H^n)} <P < 1/2 

Since 0 < a < p < 1/2 then there exists 7 > 0 such that a + p + y = 1 . From eqs. (IE. 2D one has 


PiHij)dH,j 


> 1 


fP(H,,)dHy 
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< 1 


V i,j = l,...,JV 


and 


P(Hii,Hi 2, ... ,HwN)dHiidHi2 • ■ ■ dHj, 


> 1 


P(Hii,Hi2, .. .,HNNmidH^2 ■■■dHf, 


(E.l) 


(E.2) 


(E.3) 


Now consider the systems of equations 


13 


< 1 (E. 4 ) 


and 


Y 




" f 
+ 






= 1 


= 1 


P(Hll^l2,-,tfiVN)<^^ll‘J^12-^tfiViV 


1 / 


V 




= 1 


(E. 5 ) 


(E. 6 ) 
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Eqs. HE.51) and (IE.6D represent the intersection of the unitary circle with ellipses whose major axis are equal to 






and whose minor axis are equal to 






for all i, j = 1,..., JV. Then it follows that (IE.5D and (IE.61) have solutions {u^, v^), (u, v) with u^-, v^-, u, v 9^ 0 for 
all i,j = l,...,N. Let be the eigenbasis of the Hamiltonian H. Now since a + ^ + y = 1 then one can 

define the state 


p* = a\'4>i) (V’ll+/3|'</’2) ('^2l+rl'^3> (V'sl 


(E.7) 


and the operators 


^ij = + Vq#2> ) (uij{'lp-^\+Vij{'lp2\] 




7I=(u|t/li> +v|t/) 2 ) ) (u{t/>i| +v{'l/l 2 l) (E.8) 

Since u^. +vfj = 1 , + v^ = 1 and due to the eigenbasis is orthonormal then one has that tty and it 

are projectors for all i,; = 1,... ,N. Erom Eqs. (IE.5I) and (IE.61) one obtains 

(^ij)p. = TrCp^tty) = aul + Pvf. = P{Hij)dHij V i, j = 1,... ,N 

(7i)p, = Tr(p,7r) = au^ + I3v^ = P(Hn,Hi 2 , • •. ,H^N)dH,^dH ,2 ■ ■ ■ (E.9) 

where Tr(...) denotes the trace operation. 

Due to eq. (IE.9D and since dHiidHi 2 ■ ■ ■ can be taken arbitrary small one can see that the projector 
fiij is associated with the probability of the i;-th Hamiltonian matrix element is and ft is associated with the 
joint probability of the Hamiltonian matrix elements are Hii,Hi 2 , ■ ■ ■ ,Hnn- Physically, this can be considered a 
kind of analog of the Born rule. Thus, it is reasonable to consider that there exist a relationship between ft and 
fiij for all i, j = 1,... ,N. Assuming that ft is an analytical function of the projectors one has 

00 

X! (E.IO) 

where ^kii,ki 2 ,...,kNN constant coefficients. Since the trace (n)^^ is proportional to the joint density probability 
P(Hii,Hi 2 ,• • • jHjvjv) then all the projectors Tt^ must be appear on the product in 

eq. (lE.lOD . i.e. + 0 for all i,j = 1,...,JV. Moreover, using that tty is a projector for all i = 1,... ,Af then the 
only power that survive in (lE.lOD is = 1 for all i, j = 1,... ,Af, thus one can recast (IE.101) as 

JT =!<:fill7112-(E.ll) 

where is a constant coefficient to be determined by the condition of fc is a projector. Indeed, since ft^ = ft one 
has 


ff^(^11^12 ■ ■ ■ ^Nn)^ “ i^^ll^l2 ■ ■ ■ ^NN (E.12) 

Now by taking the classical limit ft ^ 0 the product 7riifri2 • • • tinn becomes commutative so 

(^11^12 ■ ■ ■ — '^n'^12 ■ ■ ■ '^NN ~ ■ ■ ■ ^NN (E.13) 

By replacing this in eq. (IE.121) one obtains =K and since ft cannot be the null projector it follows that K = 1. 

Therefore, 71 = 71117112 ••• ftjvN- n 
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